Phantom Dark Energy Models with Negative Kinetic Term 
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We examine phantom dark energy models derived from a scalar field with a negative kinetic 
term for which V((j>) — > oo asymptotically. All such models can be divided into three classes, 
corresponding to an equation of state parameter w$ with asymptotic behavior wj, —* — 1, w$ — > 
w < — 1, and Wj, — * — oo. We derive the conditions on the potential V(<j>) which lead to each 
of these three types of behavior. For models with w<f, — > — 1, we derive the conditions on V(4>) 
which determine whether or not such models produce a future big rip. Observational constraints 
are derived on two classes of these models: power-law potentials with V(cj>) = \<ff (with a positive 
or negative) and exponential potentials of the form V(<p) = /3e A< * . It is shown that these models 
spend more time in a state with Q m ~ 0.$ than do corresponding models with a constant value of 
w^, thus providing a more satisfactory solution to the coincidence problem. 



I. INTRODUCTION 

Observational evidence Q, indicates that roughly 
70% of the energy density in the universe is in the form 
of an exotic, negative-pressure component, dubbed dark 
energy. (See Ref. Q for a recent review). If and 

are the density and pressure, respectively, of the dark 
energy, then the dark energy can be characterized by the 
equation of state parameter u>^, defined by 



w$ =P0/p0. 



(1) 



It was first noted by Caldwell Q that observational 
data do not rule out the possibility that w<f, < — 1. 
Such "phantom" dark energy models have several pecu- 
liar properties. The density of the dark energy increases 
with increasing scale factor, and both the scale factor 
and the phantom energy density can become infinite at a 
finite t, a condition known as the "big rip" 0, IU, El ■ Fur- 
ther, it has been suggested that the finite lifetime for the 
universe which is exhibited in these models may provide 
an explanation for the apparent coincidence between the 
current values of the matter density and the dark energy 
density 0. 

The simplest way to achieve a phantom model is to 
take a scalar field Lagrangian with a negative kinetic 
term. While such models have certain well-known prob- 
lems 0, El E3, El , they nonetheless provide an interest- 
ing set of representative phantom models, and they have 
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been widely studied 

In this paper, we provide a more general analysis of 
phantom models with negative kinetic terms. In the next 
section, we reexamine the behavior of such models, and 
we generalize previous studies to determine the asymp- 
totic behavior of such models and the generic conditions 
for the existence of a future singularity. In Sec. 3, we de- 
rive observational constraints on a subset of these models 
from the supernova data. In Sec. 4, we examine the so- 
lution to the coincidence problem for these models. Our 
conclusions are summarized in Section 5. 



II. PHANTOM MODELS 

A. General Results 

We limit our discussion to a spatially-flat universe, for 
which 



H*= P - 



and 



-3H(p + p), 



(2) 



(3) 



where H = a/ a, a is the scale factor, and we take H = 
c = 8itG = 1 throughout. In equations (0 and J2J, p 
and p represent the total energy density and pressure of 
the matter, radiation, and phantom fields: 



P = Pr, 

p = p r 



r Pr 

P<t>- 



(4) 
(5) 



In a phantom model with negative kinetic term, the 
energy density and pressure of the phantom are given by 



P0 = -(l/2)0 2 + K(0), 



and 



p^ = -(l/2)0 2 -V(0), 

so that the equation of state parameter is 

_ (1/2)^ + 1/(0) 
W * (1/2)0* - 1/(0) ' 

The evolution equation for this field is 
+ 3F0-F'(0) = 0. 



(6) 
(7) 

(8) 
(9) 



where the prime denotes the derivative with respect to 
(f>. A field evolving according to equation @ rolls uphill 



2 



in the potential. If the field gets stuck in a local maxi- 
mum, we simply have = — 1 and no future singularity 
S Ujl UJM ■ We will confine our attention to models for 
which V(4>) — > oo asymptotically. (A third possibility is 
a phantom model for which V(cj>) approaches a constant, 
e.g., a potential of the form V(<f>) — Vo4>/(4> + <fio)- We 
will not examine such models here, but they are discussed 
briefly in Ref. [T^l. 

Specific phantom models with V(6) — > oo have been 
previously examined in Refs. [H, [lj, Q, [H, [H, Il7| . 
while some conditions which produce a future singularity 
have been examined by Sami and Toporensky [Tsj and by 
Faraoni [l^. It is these latter studies which we will now 
generalize. 

Sami and Toporensky considered several classes of 
models, including power-law potentials, V(<f>) cx <p a with 
a > 0, exponential potentials with V(4>) ^ eA ^i an d po- 
tentials steeper than exponential (specifically, V{4>) oc 
e x<P -j _ They found that the power-law potentials asymp- 
totically produce w$ — ► — 1, with a big rip occurring only 
for the cases with a > 4. Exponential potentials lead to 
a big rip with approaching a constant (see also Refs. 



dill), while the 



,\<t> 



potential gives a singularity with 
U>0 — > — oo. 

Consider first the general relation between the asymp- 
totic behavior of V{4>) and the asymptotic behavior of 
u>0. This relation can be derived by writing the equation 
of motion as I2C 




-3(1+^) 



i + 



1 d\n(x) 
6 dln(a) 



(10) 



where is the density of the phantom field in units of 
the critical density (note that evolves with time). In 
equation I|1U|) . x = cb 2 /2V, so that x and are related 
via 



l — Wtj, 



(11) 



Equation I jlOjl differs slightly from the corresponding 
equation in Ref. j2£j because we use a different defi- 
nition of x. Equation (|10|) is the phantom version of the 
quintessence equation of motion first derived in Ref. [2l| ; 
it differs from the latter equation only in the sign of 1 +w 
on the right-hand side. 

Now consider the asymptotic evolution of in the 
limit where the universe is phantom-dominated. We can 
distinguish three possible cases, which are determined by 
the asymptotic behavior of V'/V as V — > oo. Note first 
that if V'/V asymptotically approaches a finite constant, 
a solution exists for which w approaches a constant; in 
this case we have dln(x)/d\n(a) — > in equation (jlOJI . 
On the other hand, if V'/V becomes asymptotically in- 
finite, we see that there is a solution with w — * — oo. In 
this case, x — > 1 (the latter result was noted in the nu- 
merical simulations of Ref. 18]), so that we again have 
din (a;) /din (a) — * in equation 1|10JI . Thus, from equa- 
tion HI U| l. the relation between V'/V and is given by: 



Class 1: 



Class 2: 



Class 3: 



V 

— _►()<«• to* 



(12) 



V 

— — > constant — > wq < — 1, (13) 



r 
v 



±oo <^> — > — oo. 



(14) 



These results illuminate the behavior of the specific 
examples in Refs. |l5i Il8j . For an arbitrary power-law 
potential, V oc cf> a , we have V'/V — a/ '<f>. For positive a, 
4> — > oo and V'/V — * 0, so we have u\£ — * —1, as noted in 
Ref. 0]. On the other hand, if a < 0, as in the models 
investigated below, we have <j> — > and V'/V — > — oo, so 
that W0 — > — oo. 

The only potential that corresponds exactly to Class 
2 is the exponential: V(<f>) oc e A( ^. In this case we see 
that We/, approaches a constant; equation l|l(J|) gives = 

— 1 — A 2 /3, in agreement with the results of Ref. [l5j . 
Finally, the potential V oc e A<#, ° gives V'/V = Aa</> Q-1 . 

For A > and a > 0, the field rolls in the positive <j> 
direction, and we see that V'/V — > for a < 1 (so — > 

— 1), while V'/V — >• oo for a > 1, so — > — oo. For the 
particular case a = 2, Sami and Toporensky ^| found 
Wfj, — > — oo, which agrees with our general result. 

Faraoni ^ij showed that models in which V(4>) has a 
vertical asymptote necessarily lead to a singularity in a 
finite time. All such models fall under our Class 3 (e.g., 
negative power-law potentials) , but there are also models 
in Class 3 that do not have a vertical asymptote (e.g., 
V(0) oc e A< ^° with A > and a > 1). 

All models belonging to Class 2 and Class 3 necessarily 
lead to a future singularity, but for Class 1, we must 
distinguish models that produce a big rip from those that 
do not. In the limit where — > —1, we have (p 2 /2V — > 0, 
and we can neglect the <f> term in equation so that 
equations © and J2J) reduce to fl8| 



and 



H 2 



V'{4>) 
3H ' 

_ V(<P) 
3 



(15) 



(16) 



Combining these two equations gives 



(17) 



The condition for a big rip is that V{4>) become infinite 
a finite time. From equation i|17|) . this big rip condit' 
is equivalent to the requirement that 



m 
ion 



finite 



(18) 
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FIG. 1: Evolution of the equation of state parameter Wtf, as 
a function of scale factor a (ao is the scale factor today) for 
V = \(f> 2 with fi m o = 0.25 (dot-dashed curve), fi m o = 0.3 
(dashed curve) and Q m o = 0.35 (solid curve), with (f>o = 1. 



FIG. 2: Evolution of the equation of state parameter as 
a function of scale factor a (ao is the scale factor today) for 
V = X4> 2 with n m o = 0.3, and 4>o = 0.5 (solid curve), cf>o = 1 
(dashed curve), and <j!>o = 1.5 (dot-dashed curve). 



One limit of integration in equation (|18f) is the value of 
<f> for which V — > oo, while the other limit is fixed at 
an arbitrary constant </> c . Equation (|18|l distinguishes 
models with w<j, — > — 1 that produce a big rip from those 
that do not. 

Now consider the specific example of the positive 
power-law potentials. For V(<j>) — X(f> a , with a > 0, 
the integral in equation l|18|l becomes 

roc ia/2 

L tL^"^^ (19) 

for a^4. Clearly, this integral converges (corresponding 
to a big rip) for a > 4 and diverges (giving no big rip) 
for a < 4, in agreement with the results of Ref. ^^|. The 
case a = 4 gives an integral that diverges logarithmically, 
corresponding to no big rip. Again, this agrees with Ref. 




The future singularity that occurs for Class 1 and Class 
2 has the property that the scale factor and density both 
become infinite at a finite time, t s . (In the classification 
scheme of Nojiri, Odintsov, and Tsujikawa (2^, this is a 
Type I singularity). This need not be the case for Class 3, 
in which w$ — > — oo. As noted in Ref. [is[. the potential 
V oc produces a singularity in which the scale factor 
goes to a constant at finite t s , but the density becomes 
infinite (a Type III singularity in Ref. 22]). We find 
that the negative power models investigated in the next 
section also produce a singularity with finite scale factor 
and infinite density. 



B. Specific Cases 

1. Power-Law Potentials 

Consider first the case of power-law potentials, for 
which 

v{<j>) = xr, (20) 

and we allow for the exponent a to be either positive or 
negative. Positive power-law potentials have been previ- 
ously investigated in Refs. 0, UtI ITsj . Negative power- 
law potentials have not been previously investigated in 
general, although Hao and Li [l5( showed that such mod- 
els do not lead to tracker solutions. Note that A in equa- 
tion l|20|) is chosen to fix the value of fi^ today, so it is 
not taken to be a free parameter in these models. 

We evolve these models using equation ©. We take 
= initially, but we have verified that the evolu- 
tion is independent of this initial condition. As long 
as <C p m + Pri 4> quickly decays to zero. On the 
other hand, we do find that these models are highly sen- 
sitive to the initial value for 0, which we designate as 
4>o- This contrasts sharply with the behavior of ordi- 
nary quintessence models with power-law potentials, for 
which the evolution can be insensitive to the initial con- 
ditions 0, |2j|. The actual value of 4>o is, of course, 
a model-dependent quantity. (For quintessence models, 
some estimates of plausible values of fa, produced by 
quantum fluctuations during inflation, have been derived 
by Malquarti and Liddle [24|). In this paper, we simply 
emphasize that many of our results are very sensitive to 
the assumed value of 4>q. 

We now present two representative examples, the case 
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FIG. 3: Evolution of the equation of state parameter as 
a function of scale factor a (ao is the scale factor today) for 
V — X(j> 6 with QmO = 0.3, and (f>o = 0.5 (solid curve), <fio = 1 
(dashed curve), and <f>o = 1.5 (dot-dashed curve). 



V(<f>) — \<p 2 , which does not lead to a big rip, and the 
case V(4>) — A</> 6 , which does yield a big rip. Fig. 1 shows 
the evolution of the equation of state parameter for 
V = \(j) 2 with a fixed value of <pQ (0 O = 1, in the units 
defined in Sec. II. A.), and Q m0 = 0.25,0.30,0.35 (where 
Q m o is the present-day matter density in units of the 
critical density) . Note that equation (|lf)|> implies that for 
a fixed set of initial conditions, there is a single functional 
behavior for as a function of f^. Hence, we can simply 
calculate Wrf, as a function of Q m and define the "present" 
to be the point at which Q m evolves to the given desired 
value, f2 m o- This implies that the three curves for w(a) 
differ only in an overall multiplicative constant for a; on 
the log scale in Fig. 1, the three curves are therefore 
identical and are simply displaced horizontally from each 
other. Given this simple dependence on f2 m o, and the 
fairly narrow observational limits on Q m o, we choose to 
fix the value of £l m o to be Q m o = 0.3 thoughout the rest 
of this paper. 

In Fig. 2, we show the evolution of for the same 
power-law potential for three different values of cf>o (</>o — 
0.5, 1, 1.5). Here we see a much wider divergence in the 
evolution of these models. Fig. 3 gives the evolution of 
w<j, for the power-law potential V = \<fi e and three dif- 
ferent values for <f>Q, We note that for both power laws, 
w<j, ~ — 1 at early times, but then decreases sharply as 
soon as becomes comparable to p m . Then w^, asymp- 
totically returns to —1 at late times. This behavior is 
consistent with the results for the single power-law model 
\V{4>) = ^4> 2 ] investigated in Ref. |18|. 

Now consider negative power-law potentials. In Figs. 
4 and 5, we show the evolution of for two negative 
power-law models: V — A0 _1 and V = A<fi~ 2 . As for 
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FIG. 4: Evolution of the equation of state parameter as 
a function of scale factor a (ao is the scale factor today) for 
V — A</) -1 with Q,m0 = 0.3 and 4>o — 1 (solid curve), <j> — 1.5 
(dashed curve). 



the positive power-law models, these models begin with 
Wcj, w — 1, but, as noted in the previous section, gener- 
ically evolves toward — oo. The change in behavior from 
W0 -1 to -> -co is triggered by dominating 
the expansion. Again, we see that the evolution is quite 
sensitive to the initial value of <f>. For small values of 
00 (<Po ;$ 0.5), these models yield a singularity at the 
present, a result not favored by observations. 



2. Exponential Potentials 
Now consider a second class of models of the form 

y(0)=/3exp(A<n. (21) 

For simplicity, we will consider only the case A > 0. As 
for the power-law case, it is clear from equation l|10(l that 
the functional form of u^fi^) is independent of /3, since 
V'/V is independent of (3, so we choose (5 to give the 
desired value of fi m o- 

Using our classification scheme of Sec. II. A., we find 
that 

y = Am/)"- 1 (22) 

with the field rolling in the positive <f> direction. Hence, 
the asymptotic behavior of depends on the value of a. 
For a = 1, equation l|21l) reduces to a simple exponential, 
previously investigated in Refs. 0, 0, 0] . In this case 
Wcf, ~ — 1 initially, with evolving to the value w^, 

1 — A 2 /3 at late times [l5|. As noted in the previous 
section, this is the only scalar field model with negative 
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FIG. 5: Evolution of the equation of state parameter as 
a function of scale factor a (ao is the scale factor today) for 
V — X<f>~ 2 with Q. m o = 0.3 and (f>o = 1 (solid curve), <f> = 1.5 
(dashed curve). 




FIG. 6: Evolution of the equation of state parameter as 
a function of scale factor a (ao is the scale factor today) for 
the exponential potential V = /3e A< * , with A = 1, Q. m o = 0.3, 
and d>n = 1, for the indicated values of a. 



kinetic term that corresponds exactly to Class 2 (w,p = 
wq =/= —1 at late times). For a > 1, we see from equation 
(|22|l that V /V — > oo asymptotically, so that — > -co 
at late times. (The special case of A = l,a = 2 was 
investigated in Ref. [la])- Finally, for a < 1, we see 
that V'/V — > asymptotically, so that — > —1 at late 
times. For these models, the question of whether or not 
a big rip occurs is determined by equation {THJ. We see 
that 

For < a < 1, this integral always converges, which 
indicates that for these cases we get a big rip. 

The behavior of several representative examples of 
these models is illustrated in Figs. 6 and 7. For these 
models, we take <f>o = 1, although the case a — 1 (only) 
is insensitive to the initial conditions. For all of these 
models, the transition in w^, is triggered by the onset of 
dark energy domination, as was the case for the power- 
law potentials. 



III. SUPERNOVA CONSTRAINTS ON 
PHANTOM MODELS 

Supernova la observations 0, Q provide the best cur- 
rent constraints on the value of w^. However, these limits 
normally assume either a constant w$, or a value for 
that varies linearly with a near the present. The phan- 
tom models discussed here clearly display quite different 
behavior for the evolution of w§. Hence, it is useful to 
derive supernova constraints for these models directly. 



In this section, we investigate the goodness of fit of 
various phantom models to the corresponding observed 
luminosity distance D£ fcs coming from the Snla Gold data 
set |2|. The observations of supernovae measure essen- 
tially the apparent magnitude m, which is related to the 
luminosity distance <Il by 



m(z) = M + 5\og w D L (z) , 

where 

D L (z) = —d L (z) , 
c 

is the dimensionless luminosity distance and 

d L (z) = (1 + z)d M {z) , 
with cIm(z) being the comoving distance given by 

1 



Also, 



M=M + 5 log] 



rdZ* 

,i H(z') 
c/H \ 



1 Mpc/ 



25 



(24) 



(25) 



(26) 



(27) 



(28) 



where M is the absolute magnitude. 

The data points in these samples are given in terms of 
the distance modulus 

(i ohs (z) = m(z) - M{z) =51ogd L + 25, (29) 

where is measured in Mpc. The \ 2 is calculated from 



E 



Hobs(Zi) — A*th(^; Hq, C a ) 



(30) 




FIG. 7: As Fig. 6, for A = 2. 



where the present-day Hubble parameter, Hq, is a nui- 
sance parameter that is marginalized over, and c a are the 



model parameters, 
by 



In what follows, we study A\ given 



A 2 2 2 

Ax =X -Xlcdm 



(31) 



where Xicdm ~ 177.1 is the value for the ACDM model 
with Q m = 0.3 and Q\ = 0.7, which we take as our fidu- 
cial model, as it provides a good fit to all current cos- 
mological observations. Thus, A% 2 measures how much 
worse (or better, for A\ 2 < 0) the phantom model fits 
the data compared to the standard ACDM model. 

The value of A\ 2 for the potentials V — \<fi a is shown 
in Fig. 8 as a function of a. (In this section, we take 
d>o = 1 throughout). It is clear that large positive or 
negative values of a are excluded. At the 2a level, we 
have 
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FIG. 8: The value of Ax = X ~ Xicdm as a function of a 
for the power-law potentials V = \cj> a , with Q m o = 0.3 and 
tf>o = 1. Here xlcdm ls the value of x 2 f° r the ACDM model 
with Qmo — 0.3 and Hao = 0.7. Horizontal lines give the la, 
2a, and 3a limits. 



The value of A% 2 for the exponential potentials, V = 
(3e x< ^ is shown in Fig. 9. The likelihood here is obviously 
a function of both A and a, but it is clear that smaller 
values of a are favored. Conservative limits are a < 1.5 
for A = 1 and a < 1 for A < 2. Within these limits, 
models with A = 1 allow all three classes of behavior 
for the asymptotic value of w<p, as defined in Sec. II. A. 
However, with A > 2, we see that only a < 1 is allowed, 
corresponding to — + — 1 asymptotically. As noted, 
however, these a < 1 models do produce a big rip at late 
times. 



IV. COINCIDENCE ANALYSIS 



-1.5<a<2. (32) 

A comparison of these limits with the evolution graphs 
for W0 in Sec. II. B. shows, not surprisingly, that the best 
agreement with the observations occurs for models which 
most closely resemble a cosmological constant up to the 
present. The allowed models with negative power-law 
potentials generically resemble a cosmological constant 
at early times, deviating toward large negative values 
of til,), and a singularity at late times. In this regard, 
they resemble some variants of the Chaplygin gas model 
|25|. Positive power-law models with large a, which cor- 
respond to more negative values of at the present, are 
excluded. In particular, we must have a > 4 to produce a 
big rip, and all such models are excluded by observations. 
Note, however, that we have confined our investigation 
to the initial state </>o = 1; these conclusions will differ 
for other values of <fio. 



This section is an extension of the coincidence analysis 
of Ref. |3 to the models outlined in the previous section. 
Observational evidence suggests that the matter density 
and the dark energy density have similar values today. 
However, for a cosmological constant (for instance), we 
have p\ constant, p m oc a™ 3 , so there is only a very 
small epoch in time when these two densities are within 
an order of magnitude of each other; this constitutes the 
"coincidenceproblem." 

Schcrrer |2j argued that phantom models provide a 
natural explanation of this result, since the big rip in 
such a universe is itself triggered by the onset of phantom 
domination. Hence, the universe spends an appreciable 
fraction of its total lifetime in a "coincidental" state, with 
P<f> ~ Pm- (Extensions and elaborations of this argument 
can be found in Refs. HHHH). 

As in Ref. we define the parameter r to be the ratio 
of scalar field energy density to nonrelativistic matter 
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FIG. 9: As Fig. 8, for the exponential potentials V = /3e A0 ° 
with A = 1 and A = 2. 



density: 



Pn 



(33) 



For a universe with a future singularity, we then calculate 
the fraction of the total (finite) lifetime of the universe 
that is spent in a state for which the two energy densities 
are within a factor of ro from each other: 1/ro < r < ro- 
This fraction / can be derived by first calculating the 
total lifetime of the universe, given by 



tu = 



da 



aH(a) 



(34) 



and comparing it to the time that the universe spends in 
expanding from an initial scale factor, cti, to a final scale 
factor a-2,: 



t 2 -h 



da 



oi aH(a) 



(35) 



where we choose the scale factors to correspond to energy 
density ratios of r\ = l/r for a 1 and r 2 = r for a^- The 
coincidence fraction / is then given by 



/ = t 12 /h 



(36) 



which tells us the fraction of time that the energy densi- 
ties are within a factor of ro from each other. For phan- 
tom models with constant w l p 1 this coincidence fraction 
can be expressed as 



/ = 



r(i/2) 



2w> 



ro -(2w+l)/2v 



l/r 



-dr. 



(37) 



A larger value of / corresponds to a universe in which it 
is more natural to observe ~ Q m today, ameliorating 
the coincidence problem. 
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FIG. 10: Solid curve gives the fraction of time / that the 
universe spends in a state with 0.1 < p<t,/p m < 10, for phan- 
tom models with V(<f>) = \<j> a ; f is given as a function of 
a. Dashed curve gives / for models with constant hav- 
ing the same value of at the present as the corresponding 
power-law model. 



Note that the choice of ro is somewhat arbitrary; in 
this paper we choose ro = 10, corresponding to and 
p m lying within an order of magnitude of each other. 
Since our goal is to compare scalar field phantom models 
with constant-u^ phantom models, the actual value of ro 
is not crucial. Further, we note that the value of /, like 
the evolution of the scalar field, depends on the initial 
value for <f>Q. For dcfinitcncss, we consider only </>o = 1 
in what follows, but these results are easily extended to 
other values of (f>o . 

In Fig. 10, we show the coincidence fraction / for the 
power-law models V((f>) oc <f> a as a function of a, for the 
range of values of a that is not ruled out at the 2a level 
(cf. Fig. 8). No positive values for a are displayed, since 
only the region a < 2 is consistent with the supernova 
data, while a future singularity can occur only for a > 4, 
and only models with a future singularity are relevant 
in this coincidence analysis. Within the allowed range 
for a, the value of / can be quite large, roughly 0.7 for 
the most negative allowed values of a. This means that 
the universe spends 70% of its lifetime in a "coinciden- 
tal" state, i.e., one in which the matter and dark energy 
densities lie within an order of magnitude of each other. 

Further, we note that these scalar-field-based phan- 
tom models provide a more satisfactory solution to the 
coincidence problem than the corresponding models with 
constant w$ (such as those discussed in Ref. 0) in the 
sense that the former models yield a significantly larger 
value for / for a given current value of uu. The physical 
explanation for this lies in the way in which w$ evolves in 
the models with negative power law potentials. In these 
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FIG. 11: As Fig. 10, for exponential potential V((f>) = /3e A0 ° 
with A = 1. 



models, just as in the constant-?^ models, the onset of 
the singularity is triggered by the dominance of the dark 
energy component. However, in the scalar field models, 
the beginning of phantom dominance also triggers a rapid 
decrease in w^,, producing an even more rapid evolution 
toward the singularity. 

In Fig. 11, we show the coincidence fraction / for 
the exponential models V(4>) = (3e x ^ for A = 1. Note 
that this figure encompasses all three classes of behavior 
discussed in Sec. II. A. For a < 1, we have — > — 1 
(with a big rip), for a = 1, we have w<p — > —4/3 and for 
a > 1 we have w,p — > — oo. Again, we see that over most 
of the allowed range for a, the scalar field models yield a 
larger value for / than do the corresponding constant-w^ 
models with the same present-day value for w^. 



V. CONCLUSIONS 

Our results provide a comprehensive classification of 
phantom models with negative kinetic term for which 
V{4>) — > oo asymptotically. All such models lead to ei- 
ther — > —1, Wtf, — > wo < — 1, or — > — oo, depending 
on the asymptotic form of the potential. The first set of 
these models need not lead inevitably to a big rip sin- 
gularity; we have derived the conditions on V{<p) which 
determine when a big rip occurs. Our results agree with 
all of the specific cases that have been previously inves- 
tigated elsewhere. 

Power-law potentials of the form V = \(f> a lead to 



w<j, — > — 1 for positive power laws and — > — oo for neg- 
ative power laws. However, the supernova observations 
provide both upper and lower bounds on a. For the set 
of initial conditions on (f> considered here, the only mod- 
els which are both consistent with the observations and 
lead to a future singularity are the negative power-law 
potentials. 

Generalized exponentials of the form V = fie^ can 
lead to all three sorts of behaviors indicated above. We 
have Wtj, — > — 1 for a < 1, — > wo for a — 1 and 
w<j, — > — oo for a > 1. The supernova observations favor 
smaller values of a; for the initial conditions examined 
here, we find that a < 1.5 for A = 1 and a < 1 for A = 2. 

These observational limits must be treated with some 
caution, as they are highly sensitive to the assumed ini- 
tial conditions. For the specific potentials considered 
here (power laws and generalized exponentials) we find 
that the evolution of the scalar field is very sensitive to 
the initial value of <p (but not 0), in contrast to the 
"tracker" quintessence models. Nonetheless, these lim- 
its do provide a clue to the general forms of the most 
observationally-favored models: not surprisingly, they 
are the models which most closely resemble a cosmolog- 
ical constant. 

Finally, we have considered the solution of the coinci- 
dence problem in the context of these phantom models. 
This argument relies on the idea that a universe with a 
future singularity has a finite lifetime, so it is possible to 
calculate the fraction of that lifetime, /, for which the 
densities of the matter and dark energy lie within an or- 
der of magnitude of each other. Over the observationally- 
allowed parameter range, we find that both the power-law 
potentials and the exponential potentials yield a larger 
value for / than do the corresponding constant-w^ mod- 
els (examined in Ref. 0) with the same present-day 
value for w^. Thus, these phantom models with time- 
varying Wrf, provide a better resolution of the coincidence 
problem than do constant-w^ phantom models. This is 
particularly true for models in which u>0 decreases with 
time, as such models lead to a faster evolution toward 
the singularity. 

It is clear that phantom models arising from scalar 
fields with negative kinetic terms produce a rich set of 
behaviors. The class of such models that is consistent 
with current cosmological observations yields a variety 
of different possible future fates for the universe. 
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